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Abstract 

A multicanonical update relation for calculation of the micro canonical entropy Smicro{E) by 
means of the estimates of the inverse statistical temperature /3s, is proposed. This inverse 
temperature is obtained from the recently proposed statistical temperature vi^eighted histogram 
analysis method (ST- WHAM). The performance of ST- WHAM concerning the computation of 
Smicro{E) from canonical measures, in a model with strong free-energy barriers, is also discussed 
on the basis of comparison with the multicanonical simulation estimates. 
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Reweighting and histogram methods [T]-[TO] have become essential tools for calculation 
of thermodynamic averages from samplings obtained with different coupling parameters. 
The multiple histogram method [S],[TT], also known as weighted histogram analysis method 
(WHAM), has greatly improved the efficiency of Monte Carlo simulations [121 US], mainly 
when it comes to dealing with data obtained from replica-exchange method (REM) at differ- 
ent temperatures [lll[l5]. The REM simulation has been particularly important to sample 
conformations in protein folding simulations. Generalized ensemble algorithms, in which 
REM, the multicanonical algorithm (MUCA) and their extensions P^H2U] are included, 
allow for the simulation to overcome the free-energy barriers frequently encountered in the 
energy landscapes. 

The WHAM equations combine data from an arbitrary number M of independent Monte 
Carlo or molecular dynamics simulations to enhance the sampling, thereby producing ther- 
mal averages as a continuous function of the coupling parameter, being the temperature 
the most often parameter. The solutions of the WHAM equations, obtained from energy 
data stored in histograms Ha{E), a = 1, ■ ■ ■ , M yield free-energy differences via an iterative 
numerical process. The success of the iterative process depends on the number of different 
histograms and their overlaps [iTj. More recently [21], an iteration-free approach to solve 
the WHAM equations in terms of intensive variables has been developed. This numerical 
approach, namely the statistical temperature weighted histogram analysis method (ST- 
WHAM), yields the inverse temperature Ps = dS/dE directly from a new form of WHAM 
equations. Thermodynamic quantities like entropy, can be evaluated upon numerical inte- 
gration of this statistical temperature. 

The formulation of the new WHAM equations now starts from a weighted average of the 
numerical estimates for the density of states Qa, 



with the normalization condition J2a fa{E) = 1. The density of states are estimated from 
the energy histograms H^, ^a{E) = Ha/Ha, where H^ = NaWa{E)/Za, as consequence 
of the WHAM formulation pTj. The quantities N^, Wa, and are the number of energy 
entries in the histograms H^, the sampling weights, and the (unknown) partition functions, 
respectively. The final weighted temperature is given by 
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where = HjEaHa, /3f = dlnHjdE, = -dlnWjdE, and = njE^U 



The ST- WHAM approach estimates the thermodynamic temperature from the first term 
in Eq. ^ only. This is because the second term, which amounts to the difference between 
ST- WHAM and WHAM, can be neglected for large samplings, Na » 1. This remark 
is what makes the ST- WHAM an iteration-free method to estimate the inverse effective 
temperature f3s- Now, entropy estimates follow from a careful integration of (3s{E) when 
one disregards the second term in Eq. ([3]), 



Microcanonical analysis and aims - Two sampling algorithms are successfully used 
to estimate the density of states Q{E): the Wang-Landau algorithm [221 EH] and MUCA 
[25] . This success is intimately related to their performances in obtaining a reasonable 
number of round-trips between two extremal energy values [26] . Thus, depending on whether 
the system presents strong free-energy barrier, or not, we may face some failure. Of course 
this is not a feature of these algorithms only, but in fact it is an overall behavior even for other 
generalized-ensemble algorithms [19]. Algorithms like MUCA facilitate the microcanonical 
analysis, which is important for characterization of the thermodynamic aspects of phase 
transition in small systems. Among the MUCA applications and microcanonical analysis 
are the studies of heteropolymer aggregations [271 - I30] . 

The microcanonical analysis contrasts with the usual data analysis obtained from sim- 
ulations in the canonical ensemble. For example, it is possible to observe thermodynamic 
features like temperature discontinuity and negative specific heat in the microcanonical 
ensemble [311436] . which appear at first-order phase transitions. A negative specific heat 
is a consequence of the nonconcave behavior of the microcanonical entropy Smicro{E) as a 
function of the energy and leads to the so-called convex intruder in Smicro{E) [37], a feature 
that is prohibited in the canonical ensemble. Systems where one finds convex intruders in 
the entropy present ensemble inequivalence [321 |38l [39] . It is noteworthy that the functions 
entropy S{E) and free-energy F{P) are related by the Legendre-Fenchel (LF) transform, 
= C[S{E)]. This relation is always true, independent of the shape of S{E), even with 
a nonconcave piece in its domain. The LF transform always produces a concave function. 
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On the other hand, if a function is not concave in its domain, then it cannot be obtained 
as the LF transform of another function [31]. Thus, if Smicro{E) has a convex intruder, it 
cannot be calculated from free energies in that energy domain and, as a consequence one 
has thermodynamically nonequivalent systems. With respect to WHAM, it is known that 
this method in essence provides free energies, and this has motivated our comparative study 
of the entropic results produced by ST- WHAM and MUCA. In this sense, the present study 
investigates the performance of the ST- WHAM in producing the microcanonical entropy 
from a set of data obtained with canonical measures. Another important aim here is to 
show that the integral in Eq. Q can indeed be replaced with the simple iteration relation 
used for updating the multicanonical parameters. To this end, data obtained from the two- 
dimensional (2D) Ising model with competitive nearest-neighbor ferromagnetic and long- 
range dipolar interactions are analyzed. This model presents strong free-energy barriers for 
the coupling 6 = 2, and it becomes stronger as the lattice size increases. Indeed, this model 
should be considered a benchmark for testing the performance of optimized MC algorithms. 

Methods - The multicanonical algorithm consists of sampling configurations with 
weight Wmu{E) ^ 1/Q{E). Therefore, the histogram H{E) oc Q{E)wmu{E) is almost 
flat. If the Boltzmann entropy S{E) = \nQ{E) is considered, it is convenient to write 
the following parameterization for the entropy S{E) = b{E)E — a{E), where a{E) and 
b{E) are the multicanonical parameters. Hence, the multicanonical weight is given by 
Wmu{E) = exp[—b{E)E + a{E)], with the parameter a{E) related to a multicanonical free 
energy and b{E) related to the inverse of the microcanonical temperature. 

The multicanonical parameters a{E) and b{E) are estimated from A^,. multicanonical sim- 
ulations. Usually, the number Nr is defined a posteriori when the multicanonical parameters 
present some convergence. The implementation of the multicanonical method requires the 
discretization of the energy. To this end, the labeled energies Em = Eq + me are defined, 
where each energy is associated with an integer number m. All the energies in the interval 
[Em, Em+i[ are in the mth energy bin of size e and contribute to the histogram Hmu{Em)- 
The constant Eq is defined as a reference energy near but below the ground-state energy. 
The method consists of updating the multicanonical parameters as follows. 
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b-{Em) = b--\Em) + [\nH::zJ{Em+l) ' \n H^^^J {Em)] / S , 
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where n = 1, - ■ ■ ,Nr and H^^{Em) = ma.x[ho, H^^{Ejn)], with < ho < 1. It is convenient 
to compute the above recurrence relations with the initial conditions aP{Em) = and small 
values for b'^{Em), if the simulation uses a hot-start initialization. As the multicanonical 
parameters are updated, the method samples configurations with lower energies. 

Here, the numerical integration in Eq. (|4]) is replaced with the iterative MUCA procedure 
described in Eq. (jsj): ^s{E) b{E). To illustrate this integration method, in a system 
displaying strong free-energy barriers, the 2D dipolar Ising model defined on a lattice with 
N = L"^ sites occupied by Ising spins cTj = ±1 is considered, 

n=-5J:^^^.+J:^^ (7) 

<i,j> i<j «i 

where S = J/g is the ratio between the ferromagnetic exchange interaction (J > 0) and 
the dipolar antiferromagnetic interaction of strength g > 0, and r^j is the distance between 
distinct pairs of lattice spins i and j. 

This dipolar Ising model has been employed to describe thermodynamic properties of 
ultrathin magnetic films [IQ] and presents three phases for the coupling 6 = 2 [H]. More re- 
cently, we have studied the phase transitions for this coupling using the Metropolis algorithm 
and the multiple histogram analysis 

Numerical results - Next, comparative results from ST- WHAM and MUCA simula- 
tions for L = 32 and 48 are presented. The histograms Ha were obtained at five temperatures 
Tq, and contain 3.4 x 10^ measurements. The study concerning the multicanonical results are 
based on A^^ = 10'^ multicanonical updates and 5 x 10^ histogram entries in each simulation, 
for both L = 32 and 48 lattice sizes, and with e = 1. 

The Fig. 1(a) shows the histograms obtained at Tq, for L = 32. The multicanonical 
and ST- WHAM results for the inverse temperature are depicted in Fig. 1(b). These two 
figures reveal a noisy behaviour for E/N ^ —1.18, because both the Metropolis histograms 
and MUCA do not contain appreciable energy measurements in that region. Entropy 
calculations use the ST- WHAM estimates for Ps{E) followed by the multicanonical update 
equation a{Em-i) = CL^Em) + [/^^(-Em-i) — Ps{.Em)\Em- This evaluation is denoted by ST- 
WHAM-MUCA here. Since the nonconcavity of the microcanonical entropy can hardly be 
seen on a plot, we exhibit in Fig. 1(c) the shifted microcanonical entropy ASmicro{E) = 
Smicro{E) — {A + Bs). The subtraction is performed for visualization of the nonconcavity of 
the entropy in relation to the linear function joining S^icTo{Ea) to S^icTo{Eb), where Ea and 
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Ef, are such that ASmicroiEa) = ASmicro{Eb) = 0. The MUCA estimates for Smicro{,E) are 
also included in Fig. 1(c), impressively, both estimates, ST-WHAM-MUCA and MUCA, 
yield the same results in the energy range where the phase transition takes place. This 
energy range corresponds to the so-called backbending or S-loop in the caloric curve. We 
also verified that the numerical integration of Eq. (|4]) via an analytical approximation [21] 
with energy bin size e = 1 gives the same results as the ones achieved by means of the 
ST-WHAM-MUCA procedure. 

Figure 2 displays our ST- WHAM and MUCA results for L = 48. Actually, this model 
has two phase transitions for the coupling 5 = 2, which are revealed for large lattice sizes 
only. Figure 2(a) shows the Metropolis histograms obtained at 5 temperatures. The 
results for the inverse temperature are depicted in Fig. 2(b). Backbendings are observed 
at E/N ^ —1.15 and E/N ~ —1.07, which signals the presence of phase transitions. They 
are called stripe-nematic and nematic-tetragonal, respectively. This figure shows a noisier 
behavior for the multicanonical inverse temperature as compared to the one obtained for 
L = 32, because tunneling events between the lower energy phase and the higher energy one 
are suppressed due to stronger free-energy barriers. As a matter of fact, this noisy behavior 
provides evidence for the needs of larger statistics in the MUCA simulation. Figures 2(c) and 
(d) show the shifted entropies at both phase transitions. Again, both methods give results 
in impressive agreement, mainly if we consider that the apparent disagreement between the 
ST- WHAM and MUCA results may be explained by the inefficient Metropolis sampling at 
low temperatures and in the critical region of the stripe-nematic phase transition. 

In summary, we have highlighted the performance of ST- WHAM in obtaining the mi- 
crocanonical entropy, even from only a few energy histograms produced with a canonical 
measure. A simple update procedure, based on the MUCA equations, can also be readily 
used to evaluate Smicro{E) from the inverse temperature estimates by employing ST- WHAM, 
which avoids possible numerical instabilities at a first-order phase transition. Therefore, this 
approach can be used to obtain the whole microcanonical entropy in an efficient way, and 
can stands out from other methods that also evaluate Smicro{E) but require the proper 
computation of transition matrices [43j. 
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FIG. 1. (a) histograms Ha, (b) ST- WHAM (solid line) and MUCA (dashed hne) estimates for 
the inverse temperature; (c) shifted microcanonical entropies estimated from multicanonical-like 
solution of Eq. ^ (solid line) and MUCA (dashed line) as a function of E/N for lattice size 
L = 32. 
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